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Abstract: A survey of some author’s concepts on the dynamic 
systems stability regions, in the general case of dynamic systems 
that depend on parameters is related in the paper. The property 
of separation of stable regions in the free parameters domain is 
assumed in the paper as an important property of the 
environment that is carry out and in the specified case of 
walking robot analyzed in the paper. The matrix that defines the 
linear dynamic system has the components of the matrix, 
assumed to be with real values, and the matrices that intervene in 
the exposure of the method are also, with real values of the 
components. We assumed that the matrices from the exposure 
have the complex values such that the real values are also taken 
into account as particular case of the complex values. This 
hypothesis assures a new method of analysis, in the complex 
domain, on the dynamic systems stability. Our theory on the 
stability control of the dynamic systems is applied here for 
specified walking robot model that depend on parameters. The 
critical position of the walking robot evolution is defined and 
analyzed on some cases of the walking robot leg, and possible 
application for robot walking up stairs is exposed. The further 
way of research is emphasized. 


Keywords: dynamic system; stable region; walking robots; 
critical position; walking up stair; kinematics analyze 


I. INTRODUCTION 


Any dynamic system can be considered as dynamic system 
in terms of its defining parameters without fixing their values 
as physical parameters (in particular mechanical parameters), 
geometrical parameters, possible economical, biological 
parameters or other. 

An important idea is that many concrete dynamic systems 
from the literature that depend on parameters have the proper- 
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-ty of separation between the stable and unstable regions in a 
specified domain of free parameters. The stable and unstable 
regions are separated in the domain of free parameters by a 
boundary. The property of separation can be defined by the 
fact that the stable and unstable regions are open sets, 
excepting the points on boundary. This separation aspect 
creates the freedom of stability control on a neighborhood of 
any stable point in the open stable region of the dynamic 
system. We discovered mathematical conditions of the 
stability regions existence for the dynamic systems using 
various results from applied mathematics domain. 

The property of separation of stability regions is an 
important property of the environment that can contribute to 
mathematical modelling of it. 

In this paper we analyze some aspects on the conditions of 
separation between stable and unstable regions in the free 
parameters domain of the linear or nonlinear dynamical 
systems, calling the QR algorithm applied on the real matrix 
that defines the linear dynamic system or on the “first 
approximation” of the nonlinear dynamic system, using 
Operations in the complex domain by the matrix “shift of 
origin”. The dependence of this matrix eigenvalues on the 
matrix components properties intervene in justification of 
separation property. The components of the real matrix A that 
define the linear dynamic system depending on parameters are 
assumed to be piecewise continuous in the free parameters. 


Il. ON TRANSPOSITION OF THE REAL SYSTEM IN COMPLEX ONE 
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The linear dynamic system or the “first approximation” of 
the nonlinear dynamic system depending on some parameters 
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is defined by one real matrix for the cases from the literature. 
The components of such real matrix are assumed to be 
piecewise continuous functions of the system parameters. 

The properties of the stability for the dynamic system with 
parameters are described by the properties of eigenvalues 
functions for the system matrix. The eigenvalues of the system 
matrix are the same as for the matrix in Hessenberg 
corresponding form. The Hessenberg form is defined by the 
condition a, =0 for 2<i<n, j<i-1. The dynamic 


system can be substituted by the dynamic system with the 
system matrix in Hessenberg form for the stability analysis of 
the system [28], [29]. 

The system matrix in Hessenberg form is denoted by A. 
The matrix A—A J, where / is real or complex value, is also 
a matrix in Hessenberg form. The value 2 is named “shift of 
origin” of the matrix. The shift of origin for the matrix is 
important because allow the transposition of the real matrix 
that define the dynamic system in the complex one for 
analysis of the dynamic system stability. 

The QR algorithm for matrix A, using “shift of origin”, is 
defined by the relations [29]: 


O,(A, -k,D) = Ry, Ap, = RQ; +k, =O,4,0;,5¢EN (1) 
In the above relations by A, is denoted initial matrix of the 
system in the Hessenberg form, k, is “shift of origin’, Q, is 


orthogonal matrix, R,is upper triangular matrix, A,, s22 is 


also in Hessenberg form and N is natural numbers set. 

The shift of origin, with the initial value A sufficient close 
to one initial matrix eigenvalue, real or complex, impose 
acceleration of the QR algorithm convergence to respective 
eigenvalue on the similar diagonal form of the matrix. This is 
an important property of QR algorithm with shift of origin. 

In conclusion, the real matrix A with distinct eigenvalues 
is similar with the matrix in Hessenberg form and QR 
algorithm with “shift of origin” can facilitate the convergence 
of the initial matrix to similar diagonal form of the matrix with 
real or complex eigenvalues on the diagonal. 

This study is performed in the hypothesis that all 
eigenvalues of the real matrix are distinct. For extension of the 
results in the case of real matrix multiple eigenvalues, we 
called to the study of Hirsch, Smale and Devaney that have 
demonstrated that the set of matrices with distinct eigenvalues 
from linear normed space L(R”) is opened and dense set in 
this space [1]. 

The above remark gives the possibilities to motivate the 
transmission of some properties from the real matrices with 
distinct eigenvalues to the real matrices with multiple 
eigenvalues that can intervene in evolution analysis of the 
linear dynamic system. 


Il]. ON SEPARATION OF THE DYNAMIC SYSTEM STABLE 
REGIONS 


Firstly we analyze the transmissibility of the piecewise 
continuity of the dynamic system matrix components 
depending of the system parameters to the eigenvalues of the 
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matrix. This analysis is simplified in the case of the system 
real matrix substituted by the close matrix in the complex 
domain using “shift of origin”. 


Theorem 1. If the components of the matrix Æ are 
piecewise continuous of the system parameters and the 
sequence of Hessenberg form matrices A,, s=1,2,... from 


QR algorithm that started with the matrix Æ is uniform 
convergent in the complex domain to the diagonal form of the 
matrix A then the eigenvalues of the matrix A are piecewise 
continuous. 


This theorem is verified using the property that the uniform 
convergence of the continuous functions imposes the 
continuity on the function limit. In the case in which the 
eigenvalue is on the boundary between stable and unstable 
regions in the free parameters domain this eigenvalue has null 
real part and the character of the point (stable or unstable) is 
unknown. Above theorem is capitalized here using the 
following property formulated for the continuous functions of 
one variable. 

Theorem 2. Let the function f :E >C, ECR where R 


and C are real respective complex domain, be a continuous 
function in the point x) € E and the function value f(x) 


with module | T (xo) that the 
ties æ <|f(xo)|< 2; &,BER are satisfied. Then there is a 


its SO inequali- 


neighborhood of the point x) € E where the function values 
carry out the same inequalities. 


Remark: Theorem 2 assures that the function f , 


continuous in the point x) € E preserve, in the neighborhood 
of x), the function module sign from xo. 


The mathematical conditions that assure the separation 
between stable and unstable regions for the linear dynamic 
system are described using and the following property: 


Theorem 3. Let the linear dynamic system be defined by 
the differential equation of the form d y/dt=A y(t), 


VO) HHO 0) 5 ASG, tes Sen, 
where the symbol T signifies the transposition of the matrix 


and the values a, are assumed constants. If the real part of all 


eigenvalues of the matrix A is strict negative then the solution 
of the differential equation is asymptotic stable in origin. If the 
real part of at least one eigenvalue of the matrix A is strict 
positive then the solution of the differential equation is 
unstable in origin. 

If the real parts of the matrix A eigenvalues are strict 
negative with exception of at least one eigenvalue that has null 
real part then the stability of the dynamic system in origin is 
unknown (possible stable or unstable). 


In the case of distinct eigenvalues of the real matrix A, 
attached to the linear dynamic system, a theorem on separation 
between stable and unstable regions in the free parameters 
domain of linear dynamic system, using the QR algorithm 
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with the shift of origin in complex domain is formulated 
below. This theorem is deduced calling the theorems 1- 3. 


Theorem 4 (Separation theorem). If the linear dynamic 
system defined by the real matrix A, in the Hessenberg form, 
has the piecewise continuous components of the matrix as 
functions of the dynamic system free parameters and the QR 
algorithm with the shift of origin in complex domain is 
convergent to the similar diagonal form corresponding to the 
matrix A and assures that the real part of the eigenvalue 
functions from the similar diagonal form are also piecewise 
continuous, then these conditions determine the separation 
between stable and unstable regions of the dynamic system in 
the free parameters domain. 


Remark: We comment on the possibility to substitute in 
practice the infinite QR algorithm by finite one and such that 
the conditions of the theorem 4 applications will be simplified. 

In the following we analyze some problems of walking 
robot evolution by synchronization of their leg evolution in 
kinetics hypothesis. Extension to dynamic evolution of 
walking robot, on particular case, appealing to above theory is 
intuitive possible and is very attractive to study. 


IV. MULTI-LEGGED WALKING ROBOT MODELLING 


The physical and mathematical walking robot model, with 
physical model consists from a platform and four similar legs 
that have a joined extremity attached to the robot platform and 
a synchronized evolution in the vertical plane (fig.1). 

Each of the leg is assumed to be compounded from two 
segments articulated between them and with this articulation, 
named and knee joint, traversing a circle arc route in a cycled 
evolution of the leg and where the base extremity of the leg 
traverse an imposed continuous route on the ellipse arc (fig.2). 
The leg base point in their evolution is assumed with constant 
speed component on horizontal direction, hypothesis that 
defines the leg evolution of the robot in our case studied. This 
physical model and corresponding mathematical model of the 
robot leg create the possibility to remark the possible 
appearance of the critical position, defined by the authors, on 
the circle arc of the knee joint assumed to have a trajectory on 
the circle arc in cyclic evolution of the walking robot leg. The 
points that specify the extremities of the continuous domain on 





Fig.1 Walking robot physical model. 
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circle arc where the knee joint is moving are defined in 
positions where the inferior segment of the leg is normally on 
the ellipse arc. The possibility to apply the identified critical 
positions of the walking robot evolution to climbing by stairs 
is mentioned in the paper. 

Physical model of the walking robot leg is described in 
fig.2. 


RE 4 
\\ ee Ae 
\ 
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+ 
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oO 
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Fig. 2. Physical model of walking robot leg. 

Articulated extremity of the leg to robot platform is 
denoted by Oç and the components of the robot leg are 
denoted O-Pand POQ of the 


respectively c, and where the point P is the knee joint of the 
leg with the point Q the leg base point. 


segments length a 


The knee joint P describe, in cyclic movement of robot 
leg, a trajectory on circle arc, and a base point Q of the leg an 
imposed elliptic trajectory on the superior ellipse arc O,0,. 

The circular trajectory is considered with the radius a and 
centre Oc while the elliptic trajectory is assumed with the 
semiaxes a, b and centre O,. The elliptic and respectively 


circular trajectories are transverse involving impose constant 
speed component of the leg base point, on horizontal direction. 
The critical point of the robot leg evolution is the point P. on 


the circular arc and corresponding Q, point on the elliptic arc 


with the change of the movement direction of the “knee joint”. 

The orthogonal system of coordinates for the physical 
model from fig.2 is identified by the following values of 
coordinates associated to the figure points: 


Oc (a,h), O;(a,8,), P(xp, Yp), 
O(X9; Yo), O4(x4,9), Op (xg,9). 
Some conditions are imposed on the model parameters: 
a>b>b,>0;at+c>h 


In critical position of leg “knee joint” evolution and in 
corresponding leg “base point” position is necessary null 
speed, important property that will be used in concrete robot 
leg mathematical modelled evolution. 
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The points P; and P; define extremities of the maximal 
domain on the circular arc where the point P can be moved 
because in this position the segments PQ% and P,Q, are 


normally on the ellipse arc. 

The mathematical model of the robot leg deduced from the 
physical model defined in fig.2 1s described by the implicit 
functions from (1): 


(x-a +(yp—-h)y -a° =0.; 
(x-a) /a° +(yo-h}) /b* =1. 


A covering domain for the variables from (1) is defined 
below: 


(1) 


x€[0, 2a], yp E[h-a, h+a), yo E[0, b+b]. (2) 
The explicit functions deduced from (1) have the form: 
= 2\1/2 
Yp =hF(2ax-x ) " ; 
(3) 


Yo =+ bla (2ax — x)? +b 
Let P(xp, yp) and Q(xo, Yọ) be points on the circle arc 
respectively on the ellipse arc that correspond for one position 
of the robot leg evolution. The condition (xp — xo) + 
+(Vp — Yo )’ -c* =0 is imposed. 
In the triangle O-PQ, the 


expressed as function of 0 , measure of angle (Oc PQ): 


segment OçQ can be 


(OcQ)’ =(OcP)’ +(PQ) —2(OcP\(PQ)cosO (4) 

The following analytical relation is true, where the variable 

Xo is considered an independent variable and is denoted 
by x: 

(a—x) +[h-b, F(b/a)(2ax- x°)" F = (5) 

=a* +c* —2ac cos 
The relation (5) can be used to evaluate the measure of the 
angle (O-PQ) as function of variable x, for x €[0, 2a]. 


The Q,(x,,0) and Qg(xg,0) points have abscises x, =a- 
~a (1- (b, /b) )" 7 and xp =a+a (1—-(b,/b)*)""*, where x,, 
Xg are strictly between 0 and 2a. 


The uniform linear evolution of the variable x between 0 
and 2a is assumed as below, where the selected constant 


speed œ and initial condition x) are considered: 
x(t)=Ot+X, (6) 
An evolution cycle of the robot leg can start from the 
point QO, , moving on the superior ellipse arc up to point Q,, 
using an evolution law on horizontal axe defined by (6), 
excepting a neighborhood of critical points. 
The covering domain of variables is mentioned below. 


x €[0,2a], yp Elh-a,h+a], yo €[0,5-5] (7) 


The explicit functions selected for the physical models 
from figs. 3 and 4, are as: 
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Op 


Fg.3. One physical model for walking robot leg. 


Yp =h-(2ax- x? )? 
Yo =b/a (2ax- x°)" ? -b 
In first concrete case of walking robot leg physical model is 
assured the below relation, where x €[X94, Xog]: 


(8) 


(a- x) +[(b/a)(2ax - x? ues — 


E” ae. (9) 
— (b +h) =a +c° —2ac cosé 
Because the knee joint speed is natural of zero value in critical 
positions where the sense of motion is changed, using also 
Cauchy conditions of continuity is necessary in neighborhood 
of critical points to assure continue values of displacement and 
speed as functions of time. 

In the following we describe particular cases of walking 
robot legs using physical models from the figs. 3 and 4. The 
walking robot’s physical model is defined with the help of the 
corresponding robot leg definition. 





Fig.4. Other physical model of walking robot leg. 
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For the physical robot leg model from fig. 3 the pivot point 
of leg Op, the mobile joint P of the leg that describes a circle 
arc in the leg’s evolution on cycle and the point Q of the 
robot leg base are positioned in our case on an ellipse arc of 
centre Op and semi axes length a and b. The orthogonal 


system of coordinate axes is analogue as for fig.2. 

Let Op(a,h), Op(a,—-b,), O4(% 4,9), Qg(xg,0) be points 
from figs. 3 and 4. 

An evolution cycle of the leg base point from fig.3 is 
defined by the succession Q,, Q, Qg, on the horizontal axis 


followed by the succession Q,, Q}, Q, on the superior ellipse 


arc. The critical points are not on the physical model from fig. 
3 but exist on the physical model from fig.4. 

These two models are different through the value of 
parameter A , in the first case the value A and in the second 


case the value h, > that imply the corresponding 
modification of the parameter b. 
The conditions a>b>b,>0;a+c>h are assumed, 


where the length of the inferior component of the robot leg 

is c. For the physical model from fig.3 is assumed c>a. 
Our critical position definition in the robot leg evolution, 

reminded here by us, is the position where the point P. on the 


circular arc and corresponding Q. point on the elliptic arc 


have the property of changing the movement direction of the 
“knee joint”. 

In the critical position of leg “knee joint” evolution and in 
corresponding leg “base point” evolution we remark the 
necessary null speed of leg “knee joint” and we mention again 
that this is an important property used in assurance of concrete 
evolution of the robot leg. 

The critical points P, and Q, , in our concrete case, define 


extremities of one maximal domain on the circular arc where 
the point P can be moved and in this position the 
segment P.Q, is normally on the ellipse arc. The sense of 


movement is changed on the circle arc in the extremities of the 
continuous domain of evolution for “knee joint”. In our 
opinion this notion on critical position can be generalized in 
the case more general of walking robot movement where the 
ellipse arc can be substituted by another continuous curve. 

The analytical method for identification of critical points 
on the circle arc and implicit the critical points on the ellipse 
arc, for concrete physical model from fig.4, is deduced by 
solving the equation that depend on xọ imposed by the 


orthogonally condition. 

The ellipse arc trajectory of walking robot leg base point is 
recommended by us for walking up stairs in kinematics or 
possible dynamics of walking robot evolution (see fig.5). The 
stair is configured in vertical plane, by broken line ABCD. 

The walking robot can be imagined as in fig.l. The 
trajectory of one leg, in cyclic climbing by stairs of the 
walking robot is assured ellipse arc EF, where the 
points E,F are selected critical points of the leg movement. 
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Fig. 5. Recommended trajectories of leg base point for walking up stairs 


In the critical points E, F , is necessary to be assured zero 


value of leg base point speed. The synchronized evolution of 
the legs is searched so that the stability (existence) of the 
walking robot to be assured. 

We remark that the property of walking robot stability 
(existence) position is also maintained in a neighborhood for 
one such position. 


V. CONCLUSIONS 


A mathematical study on the walking robots dynamic 
system stability control possibilities, the walking robot being 
considered as particular case of a dynamic (kinematics) 
system that depends on parameters, is performed. The 
mathematical conditions on the separation between stable and 
unstable regions, in the domain of the dynamic system free 
parameters, discovered by us, are emphasized, calling the 
algebraic operations in complex domain that permit new 
results in the stability theory. 

The property of separation described in the paper, 
encountered also in many other dynamic systems from the 
literature, without mathematical justification, is very important 
because it creates the possibility of the stability control of a 
dynamic system by optimization of the system evolution using 
the compatible criterion in the stability regions of the free 
parameters domain. The study for the dynamic systems can be 
performed, with conservation of the fundamental idea, and for 
kinematics of the walking robots. A defined walking robot 
having the property of separation in stable and unstable 
regions is analyzed. The property of stable region separation is 
described in our case of walking robot leg with the possibility 
to be extended on other cases. The critical position on the 
particular case of the walking robot is defined in the paper and 
a mathematical method of its identification is performed. The 
possibilities of our study application in control of walking 
robot evolution can be as example for walking up stairs that is 
underlined in the paper. We remark that our study has not 
exhausted the problem of dynamical systems stability control 
or for kinematics stability control and that however an 
interesting domain of walking robot stable evolution control 
has been opened. 
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